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II.  SUMMARY  OF  RESULTS 

This  project  has  centered  on  formulating  and  solving 
mathematical  problems  that  arise  in  the  study  of  systems  of 
conservation  laws  that  are  not  of  the  classical,  strictly 
hyperbolic  type.  Potential  applications  for  these  results  are 
found  in  models  for  three-phase  flow  in  porous  media,  for 
compressible  two-phase  flow,  and  for  flow  in  elastic  and 
elastoplastic  materials  (including  continuum  models  for  granular 
flow) .  Modelling  of  many  different  flow  processes  has  led  to 
systems  of  conservation  laws  in  which  the  classical  assumptions 
break  down  in  a  way  which  leads  to  distrust  of  the  models. 

Research  in  this  and  allied  projects  is  directed  at  extending  the 
mathematical  theory  of  conservation  laws.  The  practical  goal  of 
this  research  is  to  discover  which  models  are  well-posed,  and, 
hence,  to  enable  applied  scientists  to  discover  which  are  correct 
descriptions  of  various  observed  instabilities. 

The  classical  theory  of  conservation  laws  considers  a  system 

<x>  ft  u  +  I  -  0 

for  a  state  variable  U(x,t)  =  (u^,  ...  ,un)  and  flux  functions 
Fl'  *•*  '  Fm  '  Hvperbolicitv  means  that  the  eigenvalues,  A  ^ ,  of 

(2)  A ( £ , U)  =  £  qdF^U) 

are  real,  for  all  states  U  and  directions  We  are 

interested  in  systems  in  which  the  system  is  hyperbolic  in  a 
region  H  of  state  space,  but  a  pair  of  eigenvalues  of  (2) 
becomes  complex  as  the  state  U  crosses  a  boundary  2  and  enters 
a  region  8.  For  most  of  the  examples,  m  =  1  and  n  =  2,  and 
8  may  be  called  the  elliptic  region  (of  state  space) .  This 
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description  fits  models  for  dynamic  phase  boundaries  (such  as 
occur  in  two-phase  elasticity,  in  fluids  of  van  der  Waals  type  and 
in  austenitic-martensitic  transitions  in  shape-memory  alloys) ,  and 
several  models  for  low-speed  flow  (such  as  three-phase  pressure- 
driven  flow  in  porous  media  and  two-way  traffic  flow) .  The  so- 
called  single  pressure  model  for  two-phase  compressible, 
nonreacting  unsteady  flow  in  one  space  dimension  is  of  this  type 
with  n  =  4  or  6.  Current  research  on  mathematical  models  for 
granular  flows  suggests  that  change  of  type  which  is  related  to 
this  sort  of  model  may  govern  dynamic  instabilities. 

A  second  phenomenon  associated  with  systems  like  (1)  which 
describe  complex  flows  is  a  possible  lack  of  separation  of  wave 
speeds.  In  the  dynamics  of  ideal  gases,  there  is  a  single  speed, 
the  sound  speed,  which  governs  the  propagation  of  signals.  The 
subject  of  nonstrictly  hyperbolic  conservation  laws  was  motivated, 
initially,  by  consideration  of  elementary  systems,  such  as  an 
elastic  string  in  the  plane,  in  which  there  are  two  distinct  modes 
of  propagation  (longitudinal  and  transverse) ,  which  need  not  have 
distinct  speeds.  This  leads  to  the  possibility  of  nonlinear 
interactions,  resonance,  and  other  phenomena  which  may  change  the  • 
mathematical,  physical  and  computational  nature  of  the  flows. 

This  phenomenon  may  occur  in  any  flow  that  is  sufficiently  complex 
to  admit  more  than  one  characteristic  speed. 

During  the  period  of  the  AFOSR  grant,  progress  was  made  in 
the  mathematical  study  of  several  model  problems,  and  some 
underlying  principles  became  clear.  One  tentative  conclusion  is 
that  the  global  failure  of  separation  of  wave-speeds  may  be  more 
significant  in  determining  resonant  interference  and  catastrophic 
failure  of  the  model  than  is  a  local  failure  of  hyperbolicity . 

In  a  series  of  papers  [3,  7,  10],  we  study  a  system,  (a 
simple  model  for  elastic  shear  flow,)  in  which  the  wave  speeds, 
u  +  1  and  u  -  1,  are  not  globally  separated: 


v  '  13 

Vt  +  <3U  ‘  u>x  -  0 

It  is  strictly  hyperbolic  and  genuinely  nonlinear.  There  is  an 
approximate  solution  when  viscosity  is  added,  but  no  classical 
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limit  exists  as  the  viscosity  tends  to  zero.  We  find  that  a 
well-defined  but  singular  approximation  exists. 

By  contrast,  study  of  a  model  system  which  changes  type  [1, 
2],  motivated  originally  by  some  numerical  results  in  three-phase 
porous  medium  flow,  found  that  the  Riemann  problem,  even  for 
initial  states  that  are  in  the  elliptic  region,  is  "well-posed"  in 
a  sense  that  can  be  made  precise,  with  shock  waves  that  satisfy 
classical  admissibility  conditions.  The  model  system  analysed  was 
of  the  form 


(4) 


ut  +  (f(u)  -  v)x  =  0 
vt  +  (g(u) )x  =  0 


with  f  and  g  chosen  so  that  A^(u)  are  complex  for  u 
between  0  and  1,  and  so  that  and  are  both  increasing 

functions  of  u  for  u  <  0,  and  so  that  ^  is  decreasing  and 
X2  increasing  when  u  >  1.  In  [2],  it  was  shown  that  these 
assumptions  on  (4)  result  in  a  solution  of  the  Riemann  problem 
which  is  much  like  solutions  found  numerically  for  three-phase 
flow  in  porous  media  using  a  well-known  set  of  assumptions. 

Stone's  model,  for  the  three-phase  relative  permeabilities.  The 
particular  assumptions  on  the  direction  of  variation  of  the  wave 
speeds,  \ with  the  state  u  are  related  to  the  global  failure 
of  separation  of  wave  speeds  discussed  in  the  previous  paragraph, 
although  the  context,  and  the  implications  for  existence  of 
solutions  to  the  Riemann  problem,  are  completely  different  here. 

The  result  reported  in  [4]  gives  a  criterion  for  the 
appearance  of  the  type  of  Riemann  solution  which  characterizes 
this  model  (not  for  the  model  itself) .  Specifically,  there  is  an 
isolated  point,  U1#  on  the  boundary  curve  £,  where  the 
(unique)  right  eigenvector  is  tangent  to  £.  Such  a  point  is 
found,  for  example,  in  model  equations  for  three-phase  flow  in 
porous  media.  It  does  not  occur  in  models,  such  as  the  van  der 


Waals  model,  where  the  mathematically  unstable  region  corresponds 
to  a  physically  unstable  one;  one  notes  the  contrast  with  models 
in  which  the  apparent  mathematical  instability  does  not  have  a 
physical  explanation. 


We  considered  a  new  approach  to  the  question  of  admissibility 
for  shocks  in  systems  that  change  type.  One  widely  used  criterion 
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is  the  construction  of  viscous  profiles  for  different  forms  of 
viscosity:  that  is,  solutions  of 


(5) 


au 

at 


3F 

ax 


=  c  P 


2 

a^u 


ax 

which  are  functions  of  (x  -  st)/e  (  s  is  the  shock  speed). 

Here  P  is  the  viscosity  matrix.  Different  conditions  on  P  are 
appropriate  for  different  interpretations  of  the  results: 
positive  definite  or  semidefinite  symmetric  matrices  occur  in 
viscous  perturbations  of  fluid  dynamic  equations;  matrices  whose 
eigenvalues  have  nonnegative  real  part  (typically  not  symmetric) 
model  the  diffusion  introduced  by  a  first-order  finite-difference 
scheme.  We  obtained  the  following  local  result  (5]  using 
unfolding  theory  for  vector  field  dynamics:  for  a  state  UQ  in 
there  are  solutions  U, 


of  the  Rankine-Hugoniot 
which  are  in  6.  They 
states  near  one  end  can  be  joined  to  U-  by 


3(1^  -  UQ)  =  F(UX)  -  F(Uq), 


X,  near  B, 
relation, 

form  a  curve;  .  wn«  vun  ^ . —  jw.h^m  ww  wq 

forward  waves  (2-shocks) ;  near  the  other  end  there  are  backward 
profiles  (1-shocks);  there  is  an  interval  in  the  middle  of  the 
curve  where  no  profiles  exist.  This  qualitative  picture,  and  the. 
bifurcations  which  separate  the  intervals  on  the  curve,  are 
generic,  though  the  precise  location  of  the  points  which  separate 
the  intervals  depends  on  the  coefficients  of  P.  Using  this 
result,  one  can  explain  some  numerical  evidence  (noted  in  [2])  on 
competition  between  modes,  which  may  lead  to  unstable  oscillations 
in  some  cases. 

There  are  alternative  criteria  for  admissibility:  although  a 
convex  entropy  function  will  not  exist  for  a  system  which  changes 
type,  nonconvex  functions  which  have  some  physical  significance 
may  offer  a  characterization  of  admissible  shocks.  There  is  also 
a  comparison  with  steady  transonic  flow  in  two  space  dimensions: 
the  system 


( 6 )  fs  <F(U))  +  (G(U) )  =  0 

shares  many  of  the  mathematical  properties  of  (1)  when  £dF  + 
T]dG  can  have  both  real  (supersonic  states)  and  complex  (subsonic) 
eigenvalues,  and  shocks  between  the  elliptic  and  hyperbolic 
regions  are  well  known  to  occur.  Viscous  perturbation  of  this 
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equation  leads  to  a  different  mathematical  picture,  which  we 
studied  in  [8],  using  both  viscosity  and  entropy  criteria.  The 
correct  physical  admissibility  condition  (which  is  quite  different 
from  that  for  the  time-dependent  system)  is  obtained  using  real 
physical  viscosity  or  some  perturbation  of  this.  Vectorfield 
dynamics  again  isolates  an  admissibility  criterion  for  viscous 
perturbations.  Local  methods  can  be  used  to  distinguish  between 
the  steady  and  unsteady  cases;  in  [6]  it  is  shown  that  classical 
concepts  from  linear  hyperbolic  theory,  such  as  space-like 
surfaces  and  time-like  covectors,  provide  the  following 
characterization:  in  the  steady  case,  the  support  cone  of  a 

fundamental  solution  (of  the  linearized  problem)  expands  to  a 
half-plane  near  $,  while  for  unsteady  problems,  it  shrinks  to  a 
singular  ray.  Thus  it  may  be  expected,  on  purely  mathematical 
grounds,  that  the  type  of  shock  which  is  stable  will  differ  in  the 
two  cases. 

A  preliminary  study  of  the  possible  application  of  these 
techniques  to  problems  of  physical  interest  is  the  subject  of  [9]: 
the  subject  is  a  standard  model  of  two-phase  flow  which  evidences 
change  of  type.  A  simple  asymptotic  treatment  reduces  the  system  • 
to  a  pair  of  equations  which  are  of  the  same  form  as  a  classical 
model  for  two-way  traffic  flow  due  to  Bick  and  Newell  (an  early 
example  of  a  system  that  changes  type) . 

Much  of  the  work  on  change  of  type  performed  under  this  grant 
is  summarized  in  an  expository  paper,  now  in  preparation  [11]. 
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